Temperature dependence of spin correlation and charge dynamics 
in the stripe phase of high-T c superconductors 
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We examine the temperature dependence of the electronic states in the stripe phase of high-T c 
cuprates by using the t-J model with a potential that stabilizes vertical charge stripes. Charge 
and spin-correlation functions and optical conductivity are calculated by using finite-temperature 
Lanczos method. At zero temperature, the antiferromagnetic correlation between a spin in a charge 
stripe and that in a spin domain adjacent to the stripe is weak, since the charge stripe and the 
spin domain are almost separated. With increasing temperature, the correlation increases and then 
decreases toward high temperature. This is in contrast to other correlations that decrease mono- 
tonically. From the examination of the charge dynamics, we find that this anomalous temperature 
dependence of the correlation is the consequence of a crossover from one-dimensional electronic 
states to two-dimensional ones. 
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PACS numbers: 74.20.Mn, 71.10.Fd, 74.25. Jb, 78.20.Bh 



I. INTRODUCTION 

Charge stripes and related phenomena are now 
hot topics in the high-T c superconductor research 
field. Neutron-scattering measurements performed 
on Lai.475Ndo.4 Sro.i25Cu04 (LNSC) with the low- 
temperature tetragonal (LTT) structure have revealed 
the presence ofi-a charge order and an incommensurate 
magnetic order .EJ These charge and spin orders have been 
explained by assuming a stripe structure that consists of 
vertical charge stripes and antiphase between spin do- 
mains. The incommensurate magnetic order has been 
observed not only in LNSC but also in La2- :r Sr :c Cu04 
(LSC) with x ~ 0.12 carrier doping.Bu Around this car- 
rier concentration, an incomplete phase transition from 
the low-temperature orthorhorabic (LTO) phase to the 
LTT phase has been observed-ETEl While clear evidence 
of the charge order has not been reported, it is consid- 
ered that the electronic states in LSC are similar to those 
in the LNSC. 

Anomalous behaviors that are supposed to be related 
to the presence of the stripes have been observed in the 
angle-resolved photoemission (ARPES) spectrum as supu 
pressed spectral weight along the (0,0)-(7r,7r) directions 
and in the optical conductivity as the enhancement of in- 
tensity in the mid-infrared region.E3 These features have 
been explained by the present authors, by using the ex- 
act diagonalization calculation at zero-temperature for a 
model that includes both the strong electron correlation 
and the stripes, i.e., a t-J model with an additional P| 
tential introduced to stabilize vertical charge stripes 
The ground state of the model is characterized by the 
charge stripes along which the charge carriers can move 
coherently, but perpendicular to which charge carriers 
show only incoherent motion. This situation is consistent 
with jthe behavior of the Hall coefficient at low temper- 
atureB as well as the ARPES dataB in LNSC, both of 
which indicate the one-dimensional motion of the charge 



carriers.,— At the same time, the zero temperature cal- 
culationEil has shown that the spin correlation inside the 
spin domains is as strong as those for the Heiscnberg spin 
system, whereas the correlation between a spin in a spin 
domain and that in a charge stripe adjacent to the do- 
main is weak since the charge stripe and the spin domain 
are almost separated. 

With the increase of temperature, the one-dimensional 
electronic states in the stripe phase mentioned above 
are expected to be destroyed and evolved into the two- 
dimensional ones. Therefore it is interesting to know how 
such evolution occurs and how the physical quantities 
are affected by the evolution. In this study, we clar- 
ify the temperature dependence of the electronic state in 
the stripe phase. By using the finite-temperature La tic-. 
zos method developed by Jaklic and P. Prelovsek,EjE£l 
we calculate the optical conductivity and spin-correlation 
functions at finite temperatures in the t-J model with a 
potential that stabilizes the vertical charge stripes. We 
find that the spin correlation between a spin in the charge 
stripe and that in the spin domain increases with increas- 
ing temperature and then decreases toward higher tem- 
perature. This is in contrast to other spin correlations 
that decrease monotonically. By examining the depen- 
dence of charge dynamics on temperature, this anoma- 
lous temperature dependence of the spin correlation is 
concluded to be the consequence of a crossover from one- 
dimensional electronic states to two-dimensional ones. 

We introduce our model, i.e., a t-J model with a stripe 
potential, and show outlines of the finite temperature 
Lanczos method in Sec. II. In Sec. Ill, results and discus- 
sions on the temperature dependence of spin-and charge- 
correlation functions and the optical conductivity are 
presented. A possible method to confirm experimentally 
the evolution of the electronic states in stripe phase with 
temperature is also proposed. The summary is given in 
Sec. IV. 
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II. MODEL AND NUMERICAL METHOD 

We introduce the t-J model, which is given by 
H tJ = J J2 S i ■ S 3 ~ * E + Hx -) ' W 

where = Ci CT (l — rii- a ) is the annihilation operator 
of an electron with spin a at site i with the constraint 
of no double occupancy, Sj is the spin operator, and the 
summation {i,j} runs over the nearest-neighbor pairs. 

It is controversial whether .-the t-J model itself has the 
stripe-type ground stateOliia A possible origin of the 
appearance of stable stripe phase is due to the presene© 
of the long-range part of the Coulomb interactionc2rcHl 
and/or the coupling to lattice distortions. In LSCL4;he 
LTT fluctuation seems to help the latter mechanism.Qa In 
fact, the LTT structure makes Cu-0 bonds anisotropic, 
leading to directional distribution of carriers through the 
anisotropy of the Madelung potential at in-plane oxygen 
sites and that of the hopping amplitude between Cu and 
O. 

To model the directional charge distribution and the 
tendency toward the stripe instability as simply as pos- 
sible, we introduce a configuration-dependent stripe po- 
tential Vs(nh)S^ The magnitude of the stripe potential 
is assumed to depend on the number of holes rih in each 
column: Vs(nh) = for n/i=0 and 1 and Vs(n^) = — n^V 
for nh >2 with V > 0. Vs(n^) behaves like an attrac- 
tive potential for holes independent of distance between 
holes. In the following, we use a 4x4 cluster of the t-J 
model with two holes to simulate the underdoped system. 
Among four columns in the cluster, the stripe potential 
Vs{rih) is introduced into the second column from the 
left. Periodic and open boundary conditions are imposed 
along the directioiis,parallel and perpendicular to the col- 
umn, respectivelyH3 

We set parameters J/t=0A (t ~ 0.35 eV), and,V/i=l 
to obtain the ground state with charge stripes. Ej Here, 
we note that, if V = 0, the ground state in small clusters 
is not the stripe state but-a uniform state as discussed by 
Hellberg and ManousakisEJ and the energy difference be- 
tween the two states is of the order of Jo In order for the 
stripe phase to be in the ground state, it is necessary to 
introduce V with a magnitude of more than J. Therefore 
V/t=l is chosen to make the stripe state stable enough. 
Examining the V dependence of the spin and charge cor- 
relation functions discussed in the next section, we have 
found that their behaviors are not altered qualitatively if 
V > J. 

We employ nr^ 16 ^~ 

nite temperature Lanczos method.EJllj Outlines of this 
method are as follows. The statistical expectation value 
of an operator A is given by 
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where Z is the partition function defined as Z = 

L 

J2 (*z \e' l3H \ {I*/)} is a complete basis set, L is 
i=i 

the dimension of the Hamiltonian, and (3 — l/fegT, k# 
being the Boltzmann factor. Hereafter, fee is set to be 
1. In order to obtain eigenvectors and eigenvalues, we 
use the Lanczos procedure. We set an arbitrary vector 
as the initial vector of the Lanczos step \ip[), i.e., 

K) = I*;}, and calculate \ip l 2 ) , \tp l 3 ) , • • • , |V4f) ste P b y 
step, 
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where M is a given maximum number of Lanczos steps. 
Then, we obtain a tridiagonal matrix with diagonal ele- 
ments a\ with i — 1, . . . , M and off-diagonal ones b\ with 
i = 2, . . . , M, a\ and b\ being real. After diagonalizing 
the matrix, we obtain "eigenvalues" e l m and "eigenvec- 
tors" \4> l m ) ( m= 1, . . . ,M ). Since we stop the Lanczos 
steps at M steps, these "eigenvalues" and "eigenvectors" 
are approximate ones. After sampling some arbitrary 
vectors |\&;) ( I = L_ . . . , Lq; Lq <C L ) and repeating the 
above process Eq. (^), the expectation value is given by 
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L M , 2 

where Z = 2 eT^ tm \4> l m)\ ■ ^ we chose all 

1 = 1 m=l 

vectors of the basis set, i.e., L a = L and obtained eigen- 
vectors and eigenvalues exactly, Eq. (^) should be equiva- 
lent to Eq. (^). Jaklic and Prelovsek have shown that, by 

L 

using a random state \r) — ^2 Pri {flri is randomly 

i=i 

distributed) instead of \^>i), the expectation value ob- 
tained using Eq. (||) agrees well with the exact one even 
if Lq and M are much smaller than Lo Similar to static 
quantities Eq. (|J), dynamical quantities of the operator 
A are given by 
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where 



are approximate eigenvectors with approxi- 
mate eigenvalues e l n obtained by the Lanczos procedure 
Eq. (^) starting from the initial vector = A 

In small clusters, there is a characteristic temperature 
T* below which finite-size effects due to the smallness 
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FIG. 1. Hole correlation function Ch (r) in the t-J model 
without the stripe potential on a 4x4 cluster with two holes. 
J/t=0A and V/t=Q. Labeling configurations used for the 
calculation are shown in the inset. 



of the systems are appreciable. This temperature is ap- 
proximately proportional to an average level spacing in 
the low-energy sector £3 In our 4x4 t-J cluster with- 
out the stripe potential, the finite-size effects should ap- 
pear as anisotropic behaviors of physical quantities along 
the horizontal and vertical directions reflecting different 
boundary conditions along the two directions. In order 
to estimate T*, we calculate the hole correlation func- 
tion. Shown in Fig. |l| is the correlation function defined 
as 



C h (r) = (n^n)) 



(6) 



where f = Ri — Rj , Ri is the position vector for the site 
i, = 1 — rii is the hole-number operator at site i. In 
the figure, the hole correlations with two-lattice spacing 
labeled b and d are plotted, because the two correlations 
should be equivalent in magnitude in two-dimensional 
systems. The magnitude of the correlation in the b con- 
figuration is almost identical to that of d above T~0.3i. 
However, below T~0.3t, the two correlations show differ- 
ent temperature dependence. This is due to the differ- 
ence of the boundary conditions imposed on the cluster. 
Therefore the characteristic temperature T* in the t-J 
cluster without the stripe potential is about 0.3t. Even 
in the presence of the stripe potential, we assume that T* 
does not change so much from 0.3£. In fact, the average 
level spacing in the low-energy sector is similar. In the 
following, we will show the results for T/i>0.3. 

At T = 0, we employ the standard Lanczos technique 
to calculate various quantities. In this case, it is neces- 
sary to examine the dependence of given quantities on 
size and boundary condition in order to estimate finite- 
size effects. We have examined the dependence by using 
other clusters with two holes (5x4 with the same bound- 
ary condition and 4x4 with periodic boundary condition 
along both the directions) and have found that the hole 
and spin-correlation functions discussed in the next sec- 



FIG. 2. Hole correlation function d (r) in the t-J model 
with the stripe potential on a 4x4 cluster with two holes. 
J/t=0A and V/t—1.0. Labeling configurations used for the 
calculation are shown in the inset. The shadow in the inset 
denotes the direction of the vertical charge stripe. All dotted 
lines are guide to eye. 



tion depend only weakly on cluster size and boundary 
condition. 



III. RESULTS AND DISCUSSION 

Figure || shows the temperature dependence of the hole 
correlations given by Eq. (^) in the t-J model with the 
stripe potential of V/t=l. At T/t=0, the correlations of 
holes in the vertical directions labeled c and d are larger 
than those in the horizontal direction labeled a and b. 
This anisotropic behavior is due to inhomogeneous dis- 
tribution of holes induced by the stripe potential. With 
increasing temperature, the correlations in the c and d 
configurations decrease, while those in a and b increase. 
Therefore the anisotropic behavior of the correlations in- 
dicating the confinement of holes in the stripe becomes 
less pronounced at finite temperatures. The evolution 
of the hole confinement in the stripe is also observed by 
examining the hole occupation number as a function of 
temperature (not shown here); the hole number inside 
(outside) the stripe decreases (increases) monotonically 
with increasing temperature. 

Figure || shows the. results of the spin-correlation func- 
tion around a holejiil which is defined as 



C7(a,a') = ^E(^A«') 



(7) 



where a and a' denote two sites around a hole at site i in 
the stripe following the labeling convention of configura- 
tions shown in the inset of Fig. |3J, and is the number of 
holes (7Vh=2). At high temperatures, all of the nearest- 
neighbor spin correlations in the a, c, e, and / config- 
urations are antiferromagnetic and have the same mag- 
nitude, because all spins are disordered by the thermal 
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FIG. 3. Spin-correlation function around a hole C (a, a') 
in the t-J model with the stripe potential on a 4x4 cluster 
with two holes. J/t=0A and V/t=1.0. The open circle in the 
inset denotes a hole. 

fluctuation. The spin correlations between third near- 
est neighbors (the b and d configurations) are the same. 
With decreasing temperature, holes begin to be confined 
into the charge stripes and the magnitude of the spin 
correlations increase. The charge stripes make spin cor- 
relation anisotropic below T/t ~0.8. The temperature 
dependence of the spin correlations in the spin domain 
(the e and / configurations) differs from that of the spin 
correlations related with a spin in the charge stripe (a and 
c) . At the same time, the spin correlations in the b and d 
configurations become different. In the low-temperature 
region (T/t <0.4), the nearest-neighbor spin correlations 
inside the spin domains (e and f) show temperature de- 
pendence similar to that in the Heisenberg spin system 
(not shown here) because of the confinement of holes in 
the charge stripe. 

The correlation between a spin in the charge stripe 
and that in the spin domain (the a configuration) in 
Fig. U shows interesting temperature dependence below 
T/t ~0.4, when the lowest-temperature data are con- 
nected to the zero-temperature ones: The correlation is 
much suppressed at T/t—0, while other spin correlations 
(b to f) are enhanced. This anomalous temperature de- 
pendence of the a configuration can be understood in the 
following way: At T/t=0, holes move along the charge 
stripe coherently and large Drude weight is obtained par- 
allel to the stripe.Eil Spin correlation in the configuration 
a is thus suppressed in order to stabilize the coherent 
motion of holes. With increasing temperature, hole car- 
riers enter into the spin domain, and thermally fluctuat- 
ing spins in the stripe disturb the motion of holes along 
the stripe. These effects recover interaction between the 
stripe and the spin domain, and thus the spin correlation 
in the a configuration is enhanced with the gain of the ex- 
change energy. In other words, although the charge stripe 
and the spin domain are almost separated at T/t=0, the 
interaction between the charge stripe and the spin do- 
main is recovered and two-dimensional electronic system 



FIG. 4. Spin-correlation function Cs if) in the t-J model 
with the stripe potential on a 4x4 cluster with two holes. 
J/t=0A and V/t=1.0. 

is restored with increasing temperature. This picture will 
be confirmed by examining the optical conductivity as 
shown below. 

In Fig. ||, we show the result of the spin-correlation 
function given by 

Cs(r) = ^^(P(i,i)5f^) , (8) 

where r — Ri — Rj, and P (i, j) is 1 when i and j are 
in the same sublattice and it is —1 otherwise. The sum- 
mation runs over all the pairs satisfying a given r, 
and N p is the number of the pairs. As seen in the figure, 
the stripe potential, which confines holes in the charge 
stripe, makes the spin correlation anisotropic. For exam- 
ple, the spin correlation parallel to the stripe labeled c 
(d) is larger than that perpendicular to the stripe labeled 
a (6) below T/t ~0.8. With increasing temperature, the 
anisotropy decreases. 

Next, we examine the temperature dependence of the 
optical conductivity. The regular part of the optical con- 
ductivity er^f (uj > 0) ([i=x or y) is given by 

< M S M = n l ~ e ~ 0w n^ ( W ) (9) 

with 
and 

^ = (4+ M Ci -c\ci+^ , (11) 

i 

where is the eigenstate with the eigenvalue e n . e 
is the unit of the electric charge and ao is the lattice 
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FIG. 5. Temperature dependence of the regular part of the 
optical conductivity in the t- J model with the stripe potential 
on a 4x4 cluster with two holes. J/t—QA and V/t=1.0. (a) 
Perpendicular and (b) parallel to the charge stripe. Small 
broadening of 0.2t is used for delta functions. 

spacing. Other definitions are standard. Adding singu- 
lar contribution to <j r °/= (w), the real part of the optical 
conductivity reads 

(a,) = 2ne 2 D^{oj) + (w) , (12) 

where the so-called Drude weight D M is given by 



D, 



" 2N Tre 2 w 



with 



/oo 
-o 



(14) 



Here, (K^) is the kinetic energy in the /i direction. We set 
h—e—a =l hereafter. SI^uj) is calculated using Eq. (§). 

Figures ||(a) and (b) show the regular part of the opti- 
cal conductivity perpendicular, <j™ e (uj) [= a T x f(uj)], and 
parallel, ct™ s (uj) [= C7y° s (w)], to the stripe, respectively. 
At T/t=0, the intensity of a™ s (uj) is larger than that of 
a T ^ s (uj) for uj/t <1.0. In contrast, in the temperature re- 
gion of T/t > 0.3, the intensity of <r™ s (w) for u>/t <1.0 is 



0.06 

H 

Q 0.03 
0.00 > 



0.0 



\(a) —©—Perpendicular 
— «— Parallel 




— h a— Perpendicular 

— c a— Parallel 



o.o 



0.0 0.5 1.0 1.5 

Temperature T/t 

FIG. 6. (a) Temperature dependence of the Drude weight 
calculated by Eq. (jl3|) in the t-J model with the stripe 
potential on a 4x4 cluster with two holes. J/t=0A and 
V/t=1.0. (b) Temperature dependence of the integrated spec- 
tral weights I± and 1/ / of the regular part of the optical con- 
ductivity shown in Fig. U (left vertical axis) and the kinetic 
energies (K±) and (K/A (right vertical axis). In the text, 
the temperature at which I// shows a maximum is denoted 
by T m . 

much larger than that of cr™ s (w). The dramatic increase 
of a™ s (uj) from T/t=0 to T/t=0.3 is caused by the spec- 
tral weight transfer from the Drude weight parallel to the 
stripe Dj / as discussed below. Another interesting point 
in Fig. |[(b) is that the intensity of cr^ s [uS) increases with 
increasing temperature from T/i=0.3 to T/i=0.4, and 
then decreases with further increase of T. The intensity 
maximum at T/t ~0.4 is related to the maximum of In 
[= I yy ] shown in Fig. ||(b). Note that In is the inte- 

(13) 

grated spectral weight of the regular part of the optical 
conductivity a™ s (u) [see Eq. (fiA\)]. Here, let us consider 
the origin of this maximum. At T/t=0, since the motion 
of holes in the charge stripe is coherent, D//_}s very large 
while D± is almost zero as shown in Fig. @(a). On the 
contrary, I// is suppressed compared with I±. With in- 
creasing temperature, the motion of holes in the stripe 
becomes incoherent, because spins in the charge stripe 
fluctuate thermally as is evidenced by the spin correla- 
tions in the c and d configurations in Fig. |^. As a result, 
Dn is strongly suppressed, while I// is enhanced with 
the increase of temperature. The kinetic energy along 
the charge stripe (K i A [= (K y )\ exhibits small temper- 
ature dependence below T/t ~ 0.5. From the sum rule 
Eq. (H), this means that almost all of the intensity of 
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D/j transfers to the regular part 0/° s . With further in- 
crease of temperature above T/t ~ 0.5, {K//) gradually 
decreases. Since D// is almost zero in such high tem- 
perature, the decrease of {K / /) leads to the decrease of 
Therefore /// shows a maximum at the tempera- 
ture T m ~0.45i. We have examined the dependence of 
T m on the parameters V and J; T m /t~0A, 0.45, and 0.6 
for J/t— 0.2, 0.4, and 0.8, respectively, keeping V/t=1.0, 
and T m /t~0.4, 0.45, and 0.5 for V/t=0.5, 1.0, and 2.0, 
respectively, keeping J/t—OA. T m is found to be depen- 
dent not only on V but also strongly on J. These data 
are consistent with a picture that the spin degree of free- 
dom plays an important role in the charge dynamics in 
the stripe phase. 

Moreover, we find an interesting temperature depen- 
dence of the kinetic energy in Fig. |6|(b). At T/t~0, (K//) 
is larger than (Kj_) as expected. With increasing tem- 
perature, (K//) gradually decreases. Such a decrease is 
also seen in the t-J model without the stripe potential 
(not shown here). In contrast, (K±) at T/t=0.3 is larger 
than that at T/t=0. This anomalous increase of (i£±) 
is consistent with the picture proposed above, i.e., one- 
dimensional nature of the electronic states due to the 
charge stripe is destroyed by the thermal fluctuation of 
spins and thus two-dimensional electronic states are re- 
stored. 

Finally, we propose a possible method to compare the 
present results of the temperature dependence of a(oj) 
with experimental ones. The frequency dependence of 
effective carrier number N e ff(uj) is experimentally eval- 
uated from the optical conductivity by using a rela- 
tion that N e ff(uj) = 2mV/(ne 2 ) f Q a(uj')duj' , where m 
is the mass of. a free electron and V is the volume of 
the unit cell.L3 From the theoretical side, the effective 
carrier number in high-frequency limit o> — > oo is pro- 
portional to the average of the kinetic energies (K±) 
and (K//) because of a relation that / °° a(ui)duj — 
-ire 2 ((Kx) + (K//))/(4N). Since the dominant part of 
the temperature-induced change of the calculated <j//(u>) 
and a±(u)) is concentrated in the region of u <3i (~1 eV), 
the change of the averaged kinetic energies can be a 
measure of the change of the effective carrier number 
up to around 1 eV. In a realistic temperature region 
(T < J=0.4i), (K//) is almost temperature indepen- 
dent, while (K±) increases with increasing T because of 
the destruction of the stripe as mentioned above. Thus 
the averaged kinetic energy increases. We propose that 
N e ff(uj) at around w~l eV in the stripe phase of the 
high-T c cuprates increases with increasing temperature. 
This may be detectable from the detailed analysis of the 
temperature dependence of experimental data. 

IV. SUMMARY 

In summary, we have examined the temperature de- 
pendence of the electronic states with vertical charge 



stripes. The spin-correlation function, the optical con- 
ductivity, and the kinetic energy have been calculated by 
using the finite temperature Lanczos method. We have 
found that the motion of holes along the charge stripes, 
which is coherent at T/t=0, becomes incoherent with in- 
creasing temperature. We have also found that the spin 
correlation between a spin in the charge stripe and that 
in the spin domain is smaller than that for other spin 
pairs at T/t=0. However, the correlation increases with 
increasing temperature and then decreases toward high 
temperature. This anomalous increase of the correlation 
is a manifestation of the evolution to the two-dimensional 
electronic states. Moreover, we have found an anomalous 
change of the kinetic energy perpendicular to the stripe. 
All of these characteristic phenomena of the stripe phase 
appear in a realistic temperature region, i.e., T < J. The 
fact that characteristic temperatures such as T m defined 
in Fig. ||(b) depend not only on the stripe potential V 
but also on J suggests that the ordering of spins plays an 
important role in the confinement of holes in the charge 
stripe. The disordering of spins due to the thermal fluctu- 
ation destroys the charge stripe and makes the crossover 
from the one-dimensional electronic states to the two- 
dimensional ones. 
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